The sensitivity of hyperfine structure to nuclear radius and quark mass variation 
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To search for the temporal variation of the fundamental constants one needs to know dependence 
of atomic transition frequencies on these constants. We study the dependence of the hyperfine 
structure of atomic s-levels on nuclear radius and, via radius, on quark masses. An analytical 
formula has been derived and tested by the numerical relativistic Hartree-Fock calculations for Rb, 
Cd + , Cs, Yb + and Hg + . The results of this work allow the use of the results of past and future 
atomic clock experiments and quasar spectra measurements to put constrains on time variation of 
the quark masses. 
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I. INTRODUCTION 

Interest in the variation of the fundamental constants 
is motivated by theories unifying gravity and other in- 
teractions. Indications that the fundamental constants 
might have varied slightly from those of the distant 
past have been found in Big Bang nucleosynthesis and 
quasar absorption spectra (see, e.g. reviews [l|, H, 0]). 
Most publications report only constraints on possible 
variations of fundamental constants (see, e.g. reviews 
[3> El HI @)- Very stringent limits on the present time 
variation of the fundamental constants have been ob- 
tained in the atomic clock ex per iments (see, e.g. Refs. 

@, s n m mm m in mm HQ). The ma- 
jority of recent work has been devoted to the variation 
of the fine structure constant a. However, the hypothet- 
ical unification of all interactions implies that the vari- 
ation of the dimcnsionlcss strong interaction parameter 
X q = m q /K.QCD (where m q — (m u + md)/2 is the av- 
erage current-quark mass and Aqcd is the QCD scale) 
may be larger than the variation of a (see for example 
[2fj| . l2lj and the references therein) . In all intermediate 
calculations it is convenient to assume that the strong 
interaction scale Aqcd does not vary, so we will speak 
of the variation of masses (this means that we measure 
masses in units of Aqcd)- We will restore the explicit 
appearance of Aqcd in the final answers. 

In a previous paper [22| calculations of the sensitiv- 
ity of the nuclear radii to quark mass variation were 
performed. In the present paper we calculate the de- 
pendence of hyperfine transition frequency on nuclear 
radius (a preliminary approximate analytical result of 
this work was presented in Ref. [12]). Combining the 
results of the present work and Ref. [13] we calculate 
the dependence of the hyperfine structure on the quark 
masses and the dimensionless strong interaction param- 
eter X q = m q j Aqcd- These calculations are needed to 
use the results of very accurate atomic clock experiments 
to obtain constrains on the variation of the fundamental 
constants. 

The result of the present work is presented as a simple 
analytical formula. To test this formula and improve 



the accuracy we have performed numerical relativistic 
Hartree-Fock calculations for all atoms of experimental 
interest Rb, Cd + , Cs, Yb + and Hg + (and Tl in excited 
7s state). They happened to be atoms and ions with 
one s-wave electron above closed shells. One can use our 
analytical formula for other atoms where the hyperfine 
structure is dominated by s-wave electrons. For other 
electrons the effect of the nuclear radius variation is small 
and may be neglected. 



II. DEPENDENCE OF HYPERFINE 
TRANSITION FREQUENCY ON NUCLEAR 
RADIUS AND QUARK MASS 

It has been found in Ref. [22[ that the variation of a 
nuclear radius r n can be related to variation of quark 
mass by 



0.3- 



(1) 



Numerical factor 0.3 in fl} is approximately the same 
for all nuclei. Then the variation of the frequency of a 
hyperfine transition Wh may be presented as 



Sun 



where 



0.3/c/j r . 



Sr n /r n 



(2) 



(3) 



In this section we calculate khr using analytical and nu- 
merical approaches. 

To take into account finite nuclear size in the mag- 
netic dipole hyperfine structure (hfs) Hamiltonian we ap- 
proximate the nucleus by a uniformly magnetized sphere. 
Then the Hamiltonian has a form 



H. 



hfs 



U(r) 



-/j, ■ [n x a] U(r), 



l 

77 , 



r <r n 
r >r n 



(4) 
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Here n = r/r,a is Dirac matrix, \i is nuclear magnetic 
moment and r n is nuclear radius. 

For the analytical consideration we use a model of the 
nucleus where the nuclear charge is considered to be uni- 
formly distributed about a sphere of radius r n . Such a 
charge distribution corresponds to the potential 



V(r) 



r„ \ 2 2r£ J ' 



r < r„ 
r > r n 



(5) 



It is convenient to present the hypcrfinc frequency as 
tj/j = luq(1 — <5/i), where lu is the frequency at r n = 
and Sh describes the change of the hfs frequency due to 
finite nuclear radius r n . For the potential §5§ the elec- 
tron wave functions in the vicinity of the nucleus can be 
found analytically and we obtain the following approxi- 
mate expression for 5h (see appendix for details): 



S h « 7 MZr n /a B ) 2 ^ 



(6) 



where 7 = y/1 — Z 2 a 2 and a B is the Bohr radius. Then 
we obtain 



Sr n /r n 



(2j-l)5 h 
1-S h ' 



(7) 



To check these results with a more accurate approach 
we calculate atomic hfs constants using the relativistic 
Hartree-Fock method (see, e.g.[23[). Relativistic Hartree- 
Fock Hamiltonian for atoms with one external electron 
above closed shells can be written as 



H 



HF 



ca ■ p + (/? - l)mc 2 + V nuc (r) + V 



N-l 



(8) 



Here V nuc (r) is nuclear potential and V N_1 is the 
self-consistent Hartree-Fock potential of the closed-shell 
atomic core containing N — 1 atomic electrons. Nuclear 
potential V nuc (r) is found by numerical integration of the 
Fermi-type distribution of the nuclear charge. We assume 
the same electric and magnetic radius. The same Hamil- 
tonian ijHJ) is used for core and valence states. The hfs 
frequencies are expressed via expectation values of the 
hfs Hamiltonian ^ over wave functions of the valence 
electron calculated with the Hamiltonian ([8]). Note that 
neither core polarization nor correlation effects are im- 
portant for the relative change of the hfs frequency for 
s-wave energy levels since the polarization and correla- 
tion corrections are dominated by the matrix elements of 
the hyperfine interaction between the s-wave orbitals (the 
hyperfine matrix elements between the p,d,... orbitals 
are significantly smaller). We have tested this conclu- 
sion using the full-scale many-body calculations for the 
Cs hyperfine structure using approach developed in our 
work [24|. With the accuracy ~ 1% the polarization and 
correlation corrections do not change the relative value 
of the effect of the variation (to avoid misunderstand- 
ing we should note that the polarization and correlation 
corrections change the hyperfine structure constant by ~ 
40%). 



To find the change of the hyperfine frequency due to 
the change of nuclear radius we perform calculations for 
at least two different nuclear radii and than calculate 
derivative numerically 



6uj h _ uj h {r n + 5r n ) - uj h (r n - 8r n ) 



8r n 



28r n 



(9) 



The values of k^r-cai found from the Hartree-Fock calcu- 
lations are presented in Table [IJ Using ((7|) we can express 
5h in a form similar to ^ 



6 h « C(Zr n /a B ) 2 ^\ 



(10) 



where C is a fitting factor found from a comparison of 
the results of calculations with the formula f| 10[) . The 
values of C for some atoms are presented in Table HI Its 
variation from atom to atom is small and average value 
is 1.995. Using r„ = l.lAa fm and the total number of 
nuclcons A m 2.5Z leads to the formula 



6 h = 1.995 x (2.8 ■ 10- 5 Z 4 / 3 ) 2 t- 



(11) 



which can be used for any medium or heavy atom. Note 
that this result agrees with the purely analytical result 
eq. ^ to the accuracy of few per cent. 

Now we can use the expression ([2]) to calculate the 
sensitivity of uih to the quark mass due to the variation 
of nuclear radius (parameter kh q = 0.3 • khr)- All results 
are displayed in Table [Til Also presented in Table Hi] 
are the parameters K re i, k^ and k = k^ + ky iq . K re i 
is the sensitivity of the hyperfine structure to variation 
of a derived from the results of the atomic many-body 
calculations [25j]. Parameter fc M is the sensitivity of the 
nuclear magnetic moment to quark mass calculated in 
Ref. [26|. 

It is convenient to present the final results using the 
ratio of the hyperfine energy Eh = Titoh to the atomic unit 
of energy E a = m e e^ /ft 2 . Atomic experiments always 
measure the ratio of two atomic frequencies. The atomic 
unit of energy E a cancels out in such ratios. Following 
Ref. [201 we define the parameter V through the relation 



SV 



6(E h /E a ) 
E h /E a 



(12) 



Then one can use Table |TT] to find the dependence of the 
hyperfine transition frequencies on the variation of the 
fundamental constants using the following formula from 
Ref. IH : 



V 



,,2+Kr 



A, 



QCD 



(13) 



A number of limits on variation of V from different exper- 
iments are presented in Ref. [26|. These results give the 
best present time limits on the variation of m q / Kqcd- 
For example, for o;^( 87 Rb)/a;? l ( 133 Cs), we have 



A(Rb/Cs) 



U( 87 Rb) 
U( 133 Cs) 



-0.49 



ID a 



A 



QCD 



(14) 



TABLE I: The sensitivity of the hyperfine transition frequency to variation of the nuclear radius, analytical (khr) from eq. (|7I11[) 
and numerical (khr-cal) results. 



Atom or Ion gfob jj 1 Cd+ j 3 Cs ffYb+ |g 9 Hg+ g° 5 Tl(7s) ffFr 

k hr -0.010 -0.017 -0.024 -0.048 -0.077 -0.081 -0.111 

k h r-cai -0.0096 -0.0171 -0.0242 -0.0492 -0.0778 -0.0798 -0.1082 

C 1.9514 2.0034 2.0338 2.0500 2.0050 1.9623 1.9563 



and the result of measurements by [9j can be presented 
as a limit on the variation of X: 



1 dX(Rb/Cs) 
AT(Rb/Cs) It 



(-0.5±5.3)xl0- i 7yr. (15) 



Using a very stringent limit on the variation of a obtained 
using our calculations [29[ and measurements in Ref. , 

1^ = (-1.6 ±2.3) x 10- 17 /yr, (16) 
a at 

we may find the variation of X q = m q / Aqcd from eqs. 
(HI and (US]): 



1 dX q 
X~ q ~dT 



(0.3 ±2.5) x KT^/yr. (17) 



Note that the effect of the variation may be enhanced 
by 2-3 orders of magnitude in a number of molecules 
where the hyperfine splitting is approximately equal to 
an interval between the rotational levels [3011. 
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Appendix 

Analytical approach presented here is very similar to 
those used in Ref. [27| • To simplify all expressions we use 
atomic units h — e = m e = 1, c = l/a in the Appendix. 

We use an electron wave function in the form 



(18) 



where % K 3 are spherical spinors. For an s-orbital the 
initial terms of a power series solution of the radial wave 
functions inside the nucleus arc 



3 Z 2 



1 



where x — r/r n . Only those terms explicit in (|19| and 
([2TJ)l will be retained for further calculation. The external 
radial wavefunctions take the form of Bcsscl functions 



/ e =(7 + «) (J 27 {y)+bY 2l („)) 
-l{J2 1 -i{y) + bY 2l _ l {y)) 



9e =— (^2 7 (y) + bY 2j (y)) 
c 



(21) 
(22) 



where, k = {-l) j ~ l+1/2 (j + 1/2), 
7 = \J k 2 — a 2 Z 2 , y 



V8Z~i 



(23) 

As the binding energy of the electron is small in com- 
parison to the potential energy it has been neglected in 
the calculation of (|19H22[) . The constants a and b are 
found such that the wave functions remain continuous at 
r = r n . For small r n they may be approximated by 



r(2 7 ) 



' r(i + 2 7 ) 



i 



80 c 2 
-l 



+ -Zr„(3 7 + 7) 



(24) 



(25) 



As we shall only consider relative changes in the hyper- 
fine interaction normalization of the wave functions is 
not necessary. The first order correction to the energy 
is simply the s-wave expectation value of the hyperfine 
interaction which may be shown to be 



uh = {s\H hfs \s) = k j U(r)fgdr (26) 
In the limit of zero nuclear 



where k is a constant [2 
radius this becomes 

1 



UJQ 



rfogodr 



kZ 2 (l-2n) 
c( 7 — 4 7 3 ) 



(27) 



where /o = f e , go — 9e (r n — 0) ■ We define a relative change 
in the hyperfine interaction 8h by tOh = loq(1 — 5h)- 



4 



TABLE II: The sensitivity of the hyperfine structure to variation of a (parameter K T el) and to the quark mass/strong interaction 
scale m q / A.qcd (parameter k — fc M + khq) 



Atom or Ion glib ||W |fCs ffiYb+ ^ 9 Hg 4 

K ret 0.34 0.6 0.83 1.5 2.28 

k„ -0.016 0.125 0.009 -0.085 -0.088 

k hq -0.003 -0.005 -0.007 -0.014 -0.023 

k -0.019 0.120 0.002 -0.099 -0.111 



where Sg = g e — go and Sf = f e — fg. For small r n these Keeping only leading in Zr n terms reduces this to 
integrals result in 

5 h » 7 i(Zr n ) 2 ^ (30) 

g = a 2 7 (4 7 2 - 1) /J_ _ _13_^ _ ^ 

Zr n \35 1008 c 2 756 ™ Differentiation of oj/> with respect to r„ and rearranging 

i z 4 , i z 3 r„ , i 2 2 \ y ields 

-Z r„ 



640 c 4 240 c 2 360- , ^ (2 7 - l)fr fr„ 

(22^-i 27 (l + 7 )(l + 2 7 ) 6(4 7 2 -l) — = -^^T— (31) 

3r(l + 2 7 ) 2 ^ 3^Zr„ 1 j 
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